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Abstract: We present the experimental quantum tomography of 7- and 
8-dimensional quantum systems based on projective measurements in the 
mutually unbiased basis (MUB-QT). One of the advantages of MUB-QT is 
that it requires projections from a minimal number of bases to be performed. 
In our scheme, the higher dimensional quantum systems are encoded using 
the propagation modes of single photons, and we take advantage of the 
capabilities of amplitude- and phase-modulation of programmable spatial 
light modulators to implement the MUB-QT. 
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1. Introduction 

Quantum tomography (QT) is a process that allows the state reconstruction of physical systems. 
An arbitrary D-dimensional quantum system (qudit-D) is represented by a positive semidefinite, 
unit-trace Hermitian operator, which requires only — 1 independent real numbers for its 
specification. In the standard QT [1, 2], when D (i.e., a composite system of qubits), 
these parameters are determined by projecting the density operator onto completely factorized 
bases [2, 3]. The projective measurements used form a convenient set from the point of view 
of experimental implementation but, unfortunately, it does not allow the optimal quantum state 
reconstruction. Even when only local projective measurements are considered, there is another 
set that provides better performance for the QT [4]. However, as it was shown by Wooters and 
Fields [5], the optimal way for reducing the effects of the inherent statistical errors in the QT 
can only be achieved with measurements on mutually unbiased bases (MUBs), where any two 
vectors of different bases have the same overlap's absolute value. This has been demonstrated 
experimentally in [6] for polarization-entangled states and we refer to it as MUB-QT. In the 
MUB-QT, the projective measurements considered cannot be described in terms of factorized 
bases. 

Besides, while reconstructing qudit states, it is necessary to work with an overcomplete set 



of measurements that allows more accurate normalization for converting the recorded data to 
probabilities [8, 7]. In this case, the standard QT requires a number of measurements which 
scales too fast with the dimension of the quantum system. For a system composed of qubits, 
for example, it is necessary to perform 6^^ projections [9]. In the case of using MUB-QT, the 
number of measurements adopted is minimal [5, 10, 11]. For this system of qubits, it is 
necessary to perform only 2^(2^ + 1) projections. 

The advent of new technological fields, such as quantum communication, is bringing more 
motivations for the development of QT techniques. The use of qudit states has also been shown 
to be relevant for these new applications [12, 13] and, therefore, the development of techniques 
that allow better qudit state reconstructions is of upmost importance. 

In this work we present a technique that allows the MUB-QT of photonic qudit states. The 
qudits are encoded using the linear transverse momentum of single photons transmitted through 
a diffractive aperture [14, 15]. Since the qudit dimension is defined in terms of the number of 
available paths for photon transmission at this aperture, we refer to them as spatial qudits. Our 
technique relies on the fact that spatial qudit states can be coherently modified by considering 
the capability of amplitude- and phase-modulation of programmable spatial light modulators 
(SLM) [16]. This capability, combined with a spatial filtering, allows the projection of the initial 
state onto the MUBs vectors. As a figure of merit for the performance of our reconstruction, we 
use the fidelity (F) [17] of the obtained quantum states, which is defined as the overlap between 
the expected quantum state (| V^)) and the reconstructed one (p), F = (i//|p|i//). For the 7- and 
8-dimensional reconstructed qudit states, the fidelities obtained were greater than 90%. 

2. Setup description 

For a D-dimensional Hilbert space, whose dimension is a prime or a prime power number, there 
exist D+l MUBs. The constructive procedures for obtaining them have been explicitly given 
in [5, 10, 11]. The density operator of a qudit-D is represented in terms of the MUBs by [10] 

p=j: i:pif)nL«'-/, (1) 

a=lm=l 

where Ylif^ = \ i/,!"' ) ( V^,L"' |, and pif^ = Tr{pnl"^) is the probability for projecting p onto 

the m-th state 1 1//"') of the a MUB, with a = 1,2, ...,D + 1. The MUB-QT is performed by 

(a) 

determining the p;„ values from the experimental data acquired. 




Fig. 1. Experimental Setup. See the main text for details. 



The experimental setup used for the generation and reconstruction of the spatial qudit states 
is shown, schematically, in Fig. 1 . A pigtailed single mode diode laser, operating at 670 nm, is 
extremely attenuated to the single photon level by using absorptive neutral density filters (NDF). 
After this, the photons are sent to a transmissive spatial light modulator, which is composed of 
two polarizers (P) and a twisted nematic liquid crystal display (LCDl) working for amplitude- 
only modulation. This spatial light modulator is addressed with a D-multi slit as shown in Fig. 1 . 
The slit's width is 2a ~ 104 /xm and the distance between two consecutive slits \sd = 208 /im. 
The state of the single photon transmitted through the SLM can be written as [14, 16] 

i^>= E m), (2) 

1=-Id 

where /d = (D — l)/2 and the state | / ) is a single photon state defined, up to a global phase 
factor, as 

|Z) = ,P fdqe-''''''smc{qa)\lq). (3) 

\ n J 

It represents the state of the photon transmitted by the /th-slit of the SLM. 1 1^) is the Fock 
state of a photon with the transverse wave vector q. The |/) states form the logical basis in 
the D-dimensional Hilbert space of the transmitted photons, whose dimension D is defined 
by the number of slits addressed in the SLM. In this work we considered the generation of 
7- and 8-dimensional spatial qudit states. Higher dimensional qudits can, in principle, also be 
generated by considering more complexes apertures. The j3/ coefficients are dependent on the 
spatial profile and the wavefront curvature of the laser beam at the SLM plane. The first SLM 
can also be used to control each slit transmission independently, such that the initial state of 
Eq. (2) can be modified to 

\^^)„,,,, = — £ A,i3,|/), (4) 

where A/ = ^JTi, with f; representing the transmission of the /-slit, and is the normalization 
constant given by ^Y.^ff^f [16]. 

After being transmitted by the SLM, the modulated beam is imaged with a 7.5 cm focal 
length lens (LI) onto a reflective SLM, configured for doing phase-only modulation. This is 
done in a 4/-configuration, such that the LCD of the second modulator is at a distance of 30 cm 
from the first LCD. This SLM is used to modulate the phase of the beam, independently, at each 
point of sht image formation. At the plane of image formation, the diffracted photons can again 
be described as spatial qudit states [18]. Thus, the phase modulation being done by the second 
SLM allows the control of the imaginary parts of the qudit states. The state \ '^)mod\ can then 
be modified to 

1 

where 0/ is the phase given at the image of the Z-th slit. Therefore, the two SLMs allow for a 
full manipulation of the initial spatial qudit states. 

After the reflection at the second SLM, the beam is collimated by a 25 cm focal length lens 
L2, such that Li and L2 form a telescope with a magnification factor of 3.3. The collimated 
beam is then focused by a 1 m focal length lens (L3), and a point-like detector (APD) is used 
to record the single photons at the center of the interference pattern formed in the Fourier 
transform plane. The CCD camera, shown in Fig. 1, was used just for the initial alignment of the 
SLMs without attenuating the laser beam. The point-like detector is composed of a conventional 
avalanche photo-counting module with a slit of 20 /im in front of it. The spatial filtering being 
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Fig. 2. The initial qudit-7 state. In (a) one can see the aperture modulated at the first SLM. 
In (b) there is a comparison between the laser beam spatial profile and this aperture. 



introduced by the point-like detector is the last ingredient necessary for projecting the spatial 
qudit states onto the MUBs' vectors [19, 20]. The single count rate of detecting the photon at 
transversal position x at the focal plane of lens L3 (z,„— plane) is given by: 



C{x) 



(vac|£i n-*^,2m)lr>m„J2 

kxa 



(6) 



where we have used that E^^>{x,Zm) °= «Ar//3 and {w cic\ a i^^/f^\ I) — smc{kxa/f^)e '3 . is 
the focal length of lens L3. In our case, when the point-like detector is fixed at the center of the 
interference pattern formed, the single count rate will be 



C(0)^ 
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(7) 



By choosing the amplitude and phase modulation parameters as A/ — and 0/ — %J , 
respectively, we get that the single count rate at the center of the interference pattern is given 

feel (eel ^ 

by C,n °^ {ilfm \^) , where the m-th spatial vector of the a-th MUB is given by 



IWn ') = 



E4 



• (a) 



(8) 



Since in the MUB-QT an overcomplete set of measurements is considered, C,,, can be used 
to determine the probabiUty pm ' of projecting the initial state onto the m-th projector of 
the a-th MUB. This is done by changing the gray level of the apertures modulated in the first 
and the second LCD(s), respectively. For an arbitrary initial quantum state p, using the linear 

feel feel 

properties of quantum theory, we have that Q, <^ Tr(pn„, ). 

It is worth to mention that when D is a prime number, all the probabilities for a given MUB, 
can also be obtained from a single interference pattern formed by a specific phase modulation 
of the second SLM. In this case, all the probabilities of the MUB considered can be inferred by 
measuring the single counts along the transversal direction x. This means that with a single ex- 
perimental configuration one can get, only in this specific case, all the projective measurements 
of a given standard MUB labeled by the index a. 



3. Experimental results 



3.1. qudit-1 state preparation 

We now present the experimental results obtained for the MUB-QT of a spatial qudit-7 state. 
The qudit-7 state was generated with the laser beam well collimated after the optical fiber, such 
that the first SLM can be considered to be approximately at the waist of the beam. In this case 
the initial qudit state is supposed to have only real coefficients. The values of j3/ are determined 
by the spatial distribution of the beam and by the aperture addressed in the SLMl . The apperture 
modulated in this SLM is shown in Fig. 2(a). In Fig. 2(b) there is a comparison between the 
laser beam spatial profile and this aperture. For such configuration, we can expect that the 
spatial qudit state being generated in the experiment is described by | )expc = 0.256| —3) + 
0.3621 -2) +0.443I -1) +0.473|0) +0.439| 1 ) +0.352| 2) +0.254|3). 

3.2. qudit-1 state reconstruction 

The MUBs vectors for a seven dimensional Hilbert space are given explicitly by Eq. (1 1) in [5], 
and in our case they are obtained by considering at Eq. (8) that 

with a = l,2,...,7and m,Z = — 3, ...,3. (p^^^f is computed with mod(27r). The probabilities for 
projecting onto the logical base vectors, which form one of the MUBs for a 7-dimensional 




Fig. 3. Phases modulated at the second SLM and the corresponding recorded interference 
patterns. Figure (a) illustrates the relation between the gray level of the second LCD and the 
phase being modulated by this modulator. Some of the modulations used in the MUB-QT, 
for projecting the generated qudit-7 state onto | y/;" ), are shown in the insets of figures (b), 
(c) and (d). In (b) a = 7 and m = 1. In (c) a = 1 and m = 2, and in (d) a = 7 and m = 5. 
The corresponding recorded interference patterns (points) are compared with the expected 
ones (lines). The expected patterns are calculated from the amplitudes of |*I'7)eA77C' 

The 

integration time of these measurements was one second and the maximal single count rate 
recorded was 10000 counts/s. 




Fig. 4. MUB-QT of the generated qudit-7 state, (a) The expected probabiUties based on the 
predicted state | ^'j )expc- (b) The recorded probabilities with single counts, (c) The real and 
the imaginary parts of the reconstructed state. On the insets of (c) the parts of the expected 
density operator are shown. 

Hilbert space, can be determined by measuring the single counts at the plane of image forma- 
tion. This is done by removing the second SLM and scanning transversely the APD at this plane 
[18]. To generate and measure the projection probabilities for all the other seven dimensional 
MUBs vectors, it is sufficient to consider only the phase modulation implemented by the second 
SLM and the detection at the Fourier transform plane as mentioned above. In Fig. 3(a), one can 
see the dependency of the phase being modulated and the gray level of the LCD2. At Fig. 3 it 
is also shown the modulations used to project the state | ^-i)expc onto some of the vectors of the 
7-th MUB. The corresponding interference patterns were completely recorded just to show the 
purity of the modified states. The visibility of these patterns guarantees that the modulation is 
not introducing decoherence to the qudit states, showing that the modified states are nearly pure 
[16]. In Fig. 4(a) and Fig. 4(b) there is a comparison between the recorded spatial projection 
probabilities and the expected ones. The expected probabilities are calculated from the above 
given expression for | '^i)expc- One can clearly see that the obtained probabilities are very close 
to the expected ones. 

From the measured probabilities one can reconstruct the density operator using Eq. ( 1 ). How- 
ever, due to the inherent statistical fluctuations of the experimental results, the obtained density 
operator may not be positive semidefinite [2]. Thus, it is necessary to optimize it by using 
numerical techniques. For the states measured in our experiment we followed the "Forced- 
purity" approach, which was shown to provide sufficient optimization for the density operators 
of low-entropy high-dimensional qudits [21]. This is exactly our case and the qudit-7 obtained 
density operator is graphically compared with the expected one in Fig. 4(c). The expected 
real and imaginary parts are shown on the insets. The fidelity [17] between these states was 
F-jg = 0.96 ± 0.03. This error was calculated taking into account the Poissonian distribution for 
the single counts recorded in the experiment. The errors on the elements of the reconstructed 




Fig. 5. Generation and reconstruction of the spatial qudit-8 states. In (a) [(d)] there is a com- 
parison between the object being modulated at the first SLM and the spatial laser profile 
used for generating the first (second) qudit-8 state. In (b) and (e) one can see a compari- 
son between the recorded probabilities and the expected ones (insets) that are calculated 

from the states | 'Pg'' and | "Pg^' )e.t/;c, respectively. In (c) [(f)] the corresponding re- 
constructed density operator for the first (second) qudit-8 state generated is shown. On the 
insets of (c) and (f) one can see the expected density operators. 

state can be calculated directly from Eq. (1). Experimental obtained fidelities are usually limited 
by the imperfections of an experimental setup. We believe that in our case, these imperfections 
are mainly due to the mismatch that can occur while superposing the image of the first SLM 
onto the image modulated on the LCD of the second SLM. 

5.5. Preparing the qudit-S states 

In our work we also considered the generation and reconstruction of two distinct spatial qudit- 
8 states. These states were generated by changing the spatial distribution of the collimated 
laser beam at the plane of the first SLM. This was done with an optical beam expander 
mounted before the ND filters. The two spatial profiles generated for each state are shown 
in Fig. 5(a) and Fig. 5(d), respectively, where they are also compared with the object (the 8- 
multi slit) modulated at the first SLM. From these comparison we expected the first state gen- 
erated to be: = 0.2171 -|)+0.308|-f) +0.399|-|) +0.456|-i) +0.453|i) + 
0.393 1 1 ) +0.297 1 1 ) +0.202 1| ) , and the second one: 1 1'^^' = 0.343 1 - 1) + 0.350| - f ) + 
0.3551 -|)+0.358|-^)+0.359|^)+0.357| I )+0.354|§)+0.348|^). 

The case of a spatial qudit-8 state has many interesting features. The first one is that it can be 
thought as a composite system of three spatial qubits. This can easily be seen by considering the 
labeling of the I /) states as follows: ^ |000), | -§ ) ^ 1 001 ), | -| ) ^ |010), \ ^ 
|011), |i) ^ 1 100), 1 1) ^ 1 101), 1 1) ^ 1 110), 1 1) ^ I 111). These spatial qubits can, 
therefore, be used to simulate and study entangled states of three qubits. This is the case of 
the state | 'Pg'^ )expc which can not be factorized. The state | ^'"^^ )expc is almost a product state. 
In general, in case of having a 2^ dimensional qudit, with integer, this experimental setup 
allows us for mimicking a multipartite system of qubits. These subsystems are encoded 
on the transverse modes of the transmitted photon and, therefore, they cannot be used to test 



the quantum non-locality [22]. Nevertheless, they can be used for studying new techniques 
proposed for characterizing the quantum correlations of high-dimensional quantum systems 
[23, 24], or also for studying quantum information processes [25, 26, 27]. 

3.4. Reconstruction of the qudit-8 states 

Another interesting aspect of considering an qudit-8 is the fact that there are four distinct classes 
of MUBs for an eight-dimensional Hilbert space [11, 28]. The spatial qudit-8 can, therefore, 
be reconstructed by considering another type of projective measurements, that correspond to 
the projections onto the MUBs vectors of a new family of MUBs. To shown the flexibility of 
our setup, we considered the MUBs given in Table I of [11]. All the bases in this table requires 
one conditional gate acting on mimicked qubits. The corresponding MUB bases are obtained 
by applying a unitary transformation on the standard computational basis. These bases have 
the interesting properties of being optimal, when dealing with physical subsystems, since they 
require a lower number of controUed-NOT operations to be built. For sake of completeness, we 
have included all these unitary transformations in Appendix A. From these transformations we 
get the values of the {e/,"',^/,"'} parameters by using that e]"' = ^"'^ ~ ^rgl^l"^)' 

with a = 1,2,..., 8 and m,l ~ —7/2,— 5/2, ...,7/2. For generating these spatial projections it 
was necessary to change the transmission of the slits addressed at the first SLM. 

The probabilities recorded are shown in Fig. 5(b) and Fig. 5(e), where they are also compared 
to the expected ones (insets) calculated from the states | I'g'' )expc and | 'Pg^-' )ejcpc^ respectively. 
The corresponding reconstructed density operator for the first (second) qudit-8 state generated 
is shown in Fig. 5(c) [Fig. 5(f)]. The fidelities of the reconstructed states with the expected 
ones [Shown in the insets of Fig. 5(c) and Fig. 5(f)], after numerical optimization, were F^^^ = 
0.93 ± 0.03 and F^f = 0.91 ± 0.03, respectively. 

4. Conclusion 

In conclusion, we have created an experimental setup based on the use of SLMs which al- 
lows the generation and the reconstruction of distinct types of high-dimensional photonic states 
based on the mutually unbiased basis. Experimental results for the reconstruction of seven and 
eight dimensional systems were presented. The possibility of extending our study to even higher 
dimensional systems sounds feasible. For this purpose, the multi-slit arrays modulated onto the 
SLM(s), need to be replaced by more complexes apertures which will naturally define higher 
dimensional states. 

Furthermore, the presented experimental setup can be slightly modified for investigating 
other types of quantum state reconstruction of higher dimensional quantum systems. For in- 
stance, it can be used for the quantum tomography based on the informationally complete 
bases of Ref [29], unbiassed non-orthogonal bases [30], or also for the QT that is based on 
the equidistant states [31]. 

Therefore, this work can be seen as fundamental investigation which may allow more versa- 
tiUty for the experimental studies of high-dimensional quantum systems. 

A. Unitary transformations used for the qudit-8 MUB-QT 

We give the expressions for the unitary operations acting onto the logical basis | / ) of Eq. (2), 
which define the new family of MUBs used for the qudit-8 MUB-QT [11]. 
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Here, the index a in f/'"' denotes a specific mutually unbiased basis in the 8-dimensional 

[a] 

Hilbert space. It can be noted that in transformations with a = 1,2,5,8, all the j are equals to 

l/VS. Besides, in transformations with a = 4,6,7,8 we get that e^^"' = and 1/2 and finally 

for a = 3 that e^^"' = and 1 / \/2. These amplitude coefficients are implemented with the 
SLM]. Furthermore, the values for the phase modulation at the SLM2 [Eq. (8)], are given by 

(p^^ I = arg{U^^^ 1). These values belong to the discrete set of {0, 7r/2, n,3n/2}. Hence the LCD 
must allow for full modulation at the {0, 2n} domain, which is effectively provided by the used 
SLM as it is shown in Fig. 3(a). 
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